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1. Introduction 

Lanchester-type* differential-equation combat models are an important 
tool for analyzing many important problems of military operations research. 
In such a combat model, a so-called attrition-rate coefficient represents 
the fire effectiveness of a particular weapon-system type against a particular 
target type, i.e. the weapon-system type's effective firepower against such 
a target. Time-dependent attrition-rate coefficients are used to model 
temporal variations in firepower on the battlefield. Thus, we see that 
time-dependent attrition-rate coefficients are important (and, in fact, 
essential [4-6]) for the quantitative analysis of hypothetical combat. 

Militarily realistic computer-based Lanchester-type models of quite 
complex military systems have been developed for almost the entire spectrum 
of combat operations, from combat between battalion-sized units to theater- 
level operations. Nevertheless, a simple combat model may yield a clearer 
understanding of significant interrelationships that are difficult to per- 
ceive in a more complex model, and such insights can subsequently provide 
valuable guidance for more detailed computerized investigations. In this 
report we consider such a simplified variable-coefficient Lanchester-type 
model of combat between two homogeneous forces. 

For this variable-coefficient Lanchester-type model of combat 
between two homogeneous forces, different functional forms for the attrition- 
rate coefficients lead to different mathematical functions being involved 
in representing and computing the force-level trajectories. In a previous 
paper [5] we have discussed the plausibility of the hypothesis that except 


for the special case of a constant ratio of attrition-rate coefficients, 


x 
So-called after pioneering work of F. W. Lanchester [3]. 


1 


the solutions to such differential equations cannot be represented in terms 

of "elementary" functions of analysis. Thus, new transcendental functions 

arise in the study of combat modelled with time-dependent attrition-rate 

coefficients. In particular, we have previously introduced [5-6] so-called 

Lanchester-Clifford-Schlafli (LCS) functions for analyzing combat modelled 

with power attrition-rate coefficients with "no offset" (see Section 3 below). 
In the Appendix to this report is contained the most extensive 

set of tables currently available for the LCS functions: it contains tables 

of five-decimal-place values of the hyperbolic-like LCS functions F, (x), 

Hy _ 6%)» and T Cx) (see Section 4 below) for 25 fractional values of «a 

(see Section 6 below). The main body of this report provides the theoretical 

and modelling background for the use of these tables. In particular, we 

examine a model of a constant-speed attack on a static defensive position 

and show how associated range-dependent kill rates give rise to time-dependent 

attrition-rate coefficients with "no offset.'' Numerical computations are 

presented to illustrate the use of the LCS functions for analyzing es 

"aimed-fire'' combat. As a consequence of the availability of these tables, 


one can now study this variable-coefficient combat model almost as easily and 


thoroughly as Lanchester's classic constant-coefficient model. 


2. Variable-Coefficient Lanchester-Type Equations of Modern Warfare. 


We consider combat between two homogeneous forces modelled by the 


following variable-coefficient Lanchester-type [3] (see [4,5]) equations of 


modern warfare 


dx 


an -a(t)y with x(0) = Xo > 

(2,4) 
d . 
5 = -b(t)x with y(0) = Yo> 


where t = 0 denotes the time at which the battle begins, x(t) and y(t) 
denote the numbers of X and Y at time t, and a(t) and b(t) denote 
time-dependent Lanchester attrition-rate coefficients, which represent the 
effectiveness of each side's fire. These coefficients depend on variables 
such as force separation, tactical posture of targets, rate of target acqui- 
sition, firing rate, etc. (see [4-7] for further details). Variable attrition- 
rate coefficients are used to model temporal variations in firepower on the 
battlefield. In any analysis of combat, moreover, we should use the above 
equations (2.1) only for x and y > O and, for example, set dx/dt = 0 
when x = O, since negative force levels have no physical meaning. 

Mathematically, we assume that the attrition-rate coefficients a(t) 
and b(t) are defined, positive, and continuous for t, < t < +o with 


0 


ty < 0. We also assume that a(t) and b(t) € L(t),T) for any finite 


i 2 to. We further take a(t) and b(t) to be given in the form 


a(t) = k g(t) ; and b(t) = k h(t) ’ (2.2) 


where Kk and ke are positive constants chosen so that a(t)/b(t) = k/k, 


when g(t) = h(t). We introduce the combat-intensity parameter hy and 


the relative-fire-effectiveness parameter hp defined by 


i. = vk ky ; and dp = kK, : .3) 


From our assumptions about a(t) and b(t), it follows that, for example, 
T 
a(t) £ L(t) 5T) implies J. a(t)dt = +o. 
0 


The X force level as a function of time may be represented 


as [5,6] 


x(t) = x, {C, (O)C, (t) Sy (0)S_(t) J - She {Cy (0)S, (t) - S,(0)C, (t)}, (2.4) 


where the hyperbolic-like general Lanchester functions (GLF) Cy (t) and S(t) 


are linearly-independent solutions to the X force-level equation 


2 
d x 1. da{ dx _ 
a Fs dt ee a(t)b(t)x = 0 , (2:25) 
dt 
with initial conditions 
Cy (to) =l, Sy (to) = "0, 
(2.6) 


iH 
oO 
J 


{1/a(t))} dc, /dt(t)) {1/a(t))} dS, /dt(t9) =1/Vr_ . 


Here to denotes the largest finite time at which a(t) or b(t) ceases 
to be defined, positive, or continuous. The Y force level as a function 
of time is given by a similar expression, with Cy (t) and S(t) being 
analogously defined for the corresponding Y force-level equation. 


It is sometimes convenient to introduce the new independent variable 


t defined by 


cL 
t= f Va(s)b(s) ds . (2.7) 
t 
0 


It is readily seen that the transformation t = t(t) is well defined and 


invertible. Let us denote 1(0) as To: We observe that ty) < O implies 


that Tt) > 0. If we denote the "average intensity of combat" as 


va(t)b(t) , then 


ee t 
YVa(t)b(t) t = {(1/t) f Va(s)b(s) dsit =t-=- 1T 


: (2.8) 
0 0 
The substitution (2.7) transforms (2.5) into 
ds 1 d dx 
a CREE in R(t) } &¥ - x = 0, (2.9) 


with initial conditions 
and {1/YR(t,) 3 dx/dt(t,) eo. 


where R(t) = a(t)/b(t). 


3. Combat Modelled with Power Attrition-Rate Coefficients. 

The above equations (2.1) basically apply to "aimed-fire" combat 
when target-acquisition times do not depend on the numbers of targets available 
(see [5,6] for further details). A large class of tactical situations of 
interest can be modelled with the following general power attrition-rate 


coefficients [5-7] 


u _ ) 
a(t) = k(t +C)', and b(t) = k(t +C+t+A), (3.1) 


where A and C > 0. We will call A_ the offset parameter, since it allows 
us to model (with yu and vy > 0) battles between opposing weapon systems 
with different maximum effective ranges (see [5,6]). We will call C the 
starting parameter, since it allows us to model (again, with wu and vy > 0) 
battles that begin within the maximum effective ranges of the two opposing 
systems. We observe that for the general power attrition-rate coefficients 
(3.1) we have to = -C, and u and vy must be > -1l in order that a(t) 
and b(t) € L(t),T). 

The above nomenclature is motivated and possible applications of our 


work are indicated by considering S. Bonder's model of the constant-speed 


attack on a static defensive position (see [4-7] for further details) 


dx _ Gye 

cr a(r)y, and qc ee ap > qe (3¢2)) 
where r denotes the range between opposing forces, and a(r) and B(r) 
denote range-dependent attrition-rate coefficients. Range is related to time 
by 


r(t) =R,-vt , (353) 


where Ro denotes the opening range of battle and v > 0O denotes the 
constant attack speed. For example, let us consider the constant-speed attack 


of a homogeneous Y force against the static defensive position of a homo- 


geneous X force. Figure 1 diagrammatically portrays this situation. 


"JA-*Y = (31)4 Aq uaaTts st £(9)a SuoTjeaedas ao104 


*“Tepow yor}}e peeds—juejJsuod s,Japuog Jo weisetq ‘*{T ean3Ty 


aHIDG JO 
eBupy Buruedo l 


Oy aR AT IE Tee ee 


(eBudy) uolDuDdes ads04 


| 
| 
| 
| 
“(4)4 | 
| 
Ox (4) 9 an Ox 
7 | \ 
/ | 
; s 
\ 
(asuajyap 21j0;S eae: ) ! 
SHOVLLV 
SONAI30 A= Q3adS OLIV \ 
\ / 
x (4)A % / 


(4)x (4)D “eae 


The basic idea is that force separation, i.e. range between the opposing 
forces, changes over time, and the fire effectiveness of, for example, 
a single Y firer, denoted as a(r), depends on this force separation. 

In many cases of tactical interest, we may model the fire effective- 


ness of, for example, the Y weapon system (as a function of range) with 


a(r) = (34) 


where RY denotes the maximum effective range of the Y weapon system and 
u > 0. Here yu is used to model the range dependency of Y's attrition-rate 
coefficient (see Figure 2). We model 8(r) similarly, with corresponding 
quantities R, and vy _ being analogous to RO and u_ above. 


B 


If we use (3.3) to eliminate range r from (3.4), we obtain 


dx _ _ 
ans at) 3 
(3.5) 
Oy 2 
ae b(H)x 


where the time-dependent attrition-rate coefficients a(t) and b(t) are 
given by (3.1). It follows that the offset and starting parameters are given 


by 
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Bo Sie R -R 
A= ema and C= ( *«—*o), (3.6) 


and that 
wy an 
k= % (=) , and k, = Bo (=| . 7) 
a B 
We observe that A and C >O if and only if 2 ee Ro: By consider- 


ing (3.6) and Figure 3, the reader should have no trouble in understanding 
our terminology for A and C. 
When the offset parameter is equal to zero (i.e. A = 0), then the 


coefficients (3.1) reduce to 


a(t) = k (t+c)", and b(t) = k(t +c)” , (3.8) 


We will refer to (3.8) as power attrition-rate coefficients with "no offset." 


As we have seen above in Bonder's constant-speed attack model, these 
coefficients model, for example, combat between weapon systems with the same 
maximum effective range so that there is no "offset" in the "reaching out" 
of the weapon systems against each other in combat (again, see Figure 3). 
For these coefficients (3.8), the transformed X force-level equation (2.9) 


becomes 


2 
dx, (49-1) dx _ 5 2 (3.9) 


b] 


with initial conditions 
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{vt 
2 
Py pee eee (ut v+2) /2 
t= (; x ar .) (t + C) A (3211) 
and 


Gri 
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Let us observe that 0 <q <1 when u and vy >-l. Furthermore, q > 1/2 
if and only if dR/dt < 0, i.e. R(t) is a strictly decreasing function of 


time. 


4. Lanchester-Clifford-Schlafli (LCS) Functions. 


Consider the function F, (x) defined by the power series 


co 2k 
7 2) ae 
Fy@) = TC) ) attra ° Med) 
k=0 
For a #0, -l, -2,... the radius of convergence for ae 9, is infinite 


by the ratio test for convergence of power series [2]. Hence, F(z) is 


an entire function of the complex variable z= x + iy, with an essential 
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singularity at the point at infinity. Now consider the function H Gx) 


defined by the infinite series 


(xf 2° to) 
Hy@O = To) 1 er rtkee FD Oe) 
Observing that 
2a 
HG) = (1/0) (x/2) Figy ® » (4.3) 


we see that for a> 0 the infinite series (4.2) is uniformly convergent 
on compact subsets of the complex plane. From (4.3) one can readily deduce 


the recursive relation 


2 
| _(x/2)?_| Foo(x) . (4.4) 


(x) + laa + 1) | a+2 


FY Gd ~ es | 
We will call the functions F Gd) and Hx) Lanchester-Clifford-Schlafli 
(LCS) functions (see Note 10 on pp. 66-67 of [5]). Other properties are 
readily deduced and are given in Table Il. 

The function mes) satisfies the linear second-order ordinary 


differential equation 


)—-F =0, (4.5) 


with initial conditions 
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Table I. Properties of the LCS Functions Fo Gx) and H (x). 
a 


- dF /dx = Gla) 
2. dH, /dx = (x/2)°° TF (x) 
a3 FY (x) F)_4 (x) — HY GOH, _, Go) =1l #£Yx 


where @ is not an integer (including zero) 


A. FY, (x=0) = 1 

be H, (x=0) = 0 for a > 0 
6. dF /dx(x=0) = 0 

7. {(«/2)"*% dH /dx} ig = 1 

8. Fi s2®) = cosh x 

g, Hy 20%) = sinh x 
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dF 
a 


F (0) =1, and ie (0) =0 , 
while H Gx) satisfies 
aH, 2a-L a 
9 a ( x ) dx = Hy = 0 > (4.6) 
dx 


with initial conditions 


1-2, dH 
H(O) = 0, and e a =1., 
x=0 


Thus, {Fs Hy 4! is a fundamental system of solutions to 


cf+ (el Z_peo, (4.7) 


with Wronskian W(F Hy _4) = Gar. It follows that the GLF for the 
X and Y force-level equations for combat modelled with the attrition-rate 


coefficients (3.8) are given by 


| hy \" 
| Cy (t) = Hite): Sy (t) a ee ; (4.8) 
| 
dy 1-2q 
Cy (t) = ce De Sy (t) hes] H Retey? ; (4.9) 


where p = 1-q. If we define 


i: 


T(x) = HL GO/F Gx) (4.10) 
then 
S, (t) ( dy oe H (t(t)) 
ene CS ae Oe ae F Gr)’ Reet) 
X q 
or 
( dy ia 
T(t) = hae) oes . (4.12) 


where TY (t) denotes a hyperbolic--like GLF, which corresponds to the 


hyperbolic tangent. Observing that for u, v > -l, lim t(t) = +, 


Tre +o 


we see that T (x) is a strictly increasing function of x on the interval 


[O, + ~) and 


0<T (x) < +o. for O<x<+e, (4.13) 
with 
T (1-a) 
Lin» Tote) = P (4.14) 
ete, +2 T (a) 


since by the results of Taylor and Comstock [7] the parity-condition parameter 


Q* = Q¥(u, v, C = 0) is given by 


r 2q-1 
7 7 1 - I r(p) 
Lim Ty (t) = O*Cu, v, 0) = sea) Tq) (4315) 


t>r~ +2 


We recall that Taylor and Comstock [7] have introduced the so- 
called parity-condition parameter Q* as the value (or range of such values) 


for the initial condition Q _ to the initial-value problem 


16 


= 1 7 : - 2 
a — a(t)E, with Ey (tq) =l1l, 
. (4.16) 
dE 
XY. . - a 7 
ae Vip b(t)E, with E,(t,) =Q, 


such that Ey (t3Q*) and H.(t30*) = 0 for ail. 6 Se In other words, 


0° 
Q* is the value of Q in (4.16) above such that neither E, nor E 
ever become zero. In this case, both Ey (t3Q*) and E. (t3Q*) are positive, 
strictly decreasing functions, similar to decreasing exponentials. Thus, 
we may call Q* "the Y equivalent of an X force of unit strength," 
since the forces are "at parity," with neither force being annihilated in 
finite time. Taylor and Comstock have shown that for either a(t) f L(0,+ ~) 
or b(t) f L(O, + ~), then Q* is unique and given by 

S(t) 


fo. ae Se Ct (4.17) 
topo y(t) Oo 


The significance of the parity-condition parameter Q* is that it allows 


us to predict force annihilation as the following theorem shows. 


THEOREM 1 (Taylor and Comstock [7]): Assume that either 
a(t) f L(O, +”) or b(t) f L(O, +”). Then the X _ force 


will be annihilated in finite time if and only if 


| x C(O) - Q*S_.(0) 
0 xX xX 
Yo Vip oe = S,, (0) : oe 


L? 


5. Use of LCS Functions for Analyzing Combat. 
The Lanchester-Clifford-Schlafli (LCS) functions F(x) and H, (x) 


are useful for analyzing "aimed-fire" combat (see Section 3 above) modelled 
with the power attrition-rate coefficients with "no offset" (3.8), which 


we rewrite here as 
a(t) = k(t +c)", and b(t) = k(t +c)”. (o:1) 


In other words, the LCS functions arise in solving the differential combat 
model (2.1) with attrition-rate coefficients (5.1). In order that both 
a(t): and) b(t) € L(t),T), we must have u and v > -l. Military situations 
modelled by these equations have been discussed in Section 3 above, e.g. 
combat between two weapon systems with the same maximum effective range. 
For such combat, the LCS functions may be used to 

(1) compute force-level declines, 

(2) predict force annihilation,” 
and (3) predict the time of force annihilation. 

Let us now see how the LCS functions may be used to obtain the 
above information about force-level declines and force-annihilation pre- 
diction. According to (2.4), (4.8), and (4.9) above, the xX force level 


is given by 
x(t) = Xo {F (1) F(r(t)) - rar H Cees 


- 


2q-1 
= Fo Ve (2s) CRT a HOVER oH Wig ete) (S22) 
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where q is given by (3.10), p = l-q, and ,;(t) is given by (3.11), which 


we rewrite as 


2 


The time to annihilate the X force* is determined by (te) = 0, and 


it follows that 


d q-p 
x,F (t,) + y ne i Cae 
T Ea Cr a Ne 0 ie ° (5.4) 
qa A; = \a-P 
XH, To) + Yo flees) F{%? 
where from (4.10) 
T (7(t)) = 3 A7(b)) /F rte). (5.5) 
q Pp q 


and we recall that p+q#=t1. It follows that the time to annihilate X, 


xX , ; 
t is given by 


* 
If we multiply the first equation of (2.1) by y, the second by x, add, 
and integrate, we obtain 


t 


x(t) y(t) = Xoo 7 f {a(s) ote) + b(s) ide ; 
0 


which shows that x(t) and y(t) can have at most one finite zero. Hence, 
if x(tX) = 0, then we know that y(t) > O for all t > 0. 
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(t,) / ( a ie ( 
~ 7) 2,1 x,F. Ct) + ¥, YAS | —— > H_(1,) 
a) es a a o’"’R \p+tvt2 po (ea 


x q-P 
I 
Xgl? + Yo Vp @ rere aay 


Taylor and Comstock [7] have shown that Pe, is strictly increasing and 


satisfies (see also (4.12) above) 
cS < [(p)/T(q) , (Sa7) 


where p = 1l-q. It follows that in order for X to be annihilated in finite 


time, the right-hand side of (5.4) must be less than [(p)/I'(q). Let us 


observe that for to = -C = 0, (5.4) simplifies to 
P-q 
r 
T(t) = esi) —s : (5.8) 
¥p Op 


Thus, we have proved the following theorem concerning force- 


annihilation prediction. 


THEOREM 2: Consider combat between two homogeneous forces 
modelled by (2.1) with power attrition-rate coefficients 
(5.1). Assume that yu and vy > -1 and that the above 
equations hold for all time. Then the X _ force will be 


annihilated in finite time if and only if 
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q-P 
= | 
ine? | XFL (to) ‘J Yo Ap (. $y Z Hy) 


q-P 

nN 

| I | 

< PO) ) XoHy (tg) + YQ Mg ‘ Fu+2 FS? 3 eae 
where q = (v +1)/(u +v +2) and p = 1-q. For to = 0 

(i.e. C= 0 so that io 0), X will be annihilated in finite 


time if and only if 


Xo r(p) ( Ay \ 
_Y . sp) a ee 
Yo l'(q) Vip utvt2 , (5.10) 


6. Tabulation of LCS Functions. 

This report contains the most extensive set of tables of the 
Lanchester-Clifford-Schlafli functions currently available. The Appendix 
contains tables of five-decimal-place values of the hyperbolic-like LCS 


functions F(x), H > and T(x) for various values of the argument x, 


e. 
namely x = 0.00 (0.01) 2.00 (0.1) 10.0, and a = 1/2, 1/3, 2/3, 1/4, 3/4, 
Ny 26 9154 GI 35 2015. 317 6 SIs: BI: S795. O19¢ Slits. S/ik, G/ity. BL, 
5/13, 8/13, 5/17, 12/17, 5/21, and 16/21. These values of the index «a 


5, 2, and 3 in (3.8) and allow one 


to analyze, for example, a fairly wide variety of range capabilities for 


correspond to u, v = 0, 1/4, 1/2, 1, 1 


weapon systems in the constant-speed-attack model of Section 3. These 


waa 


tables have been calculated by the recursive means given in Section 8 of 
[5]. A less extensive tabulation (namely, for a = 1/2, 1/3, 2/3, 1/4, 3/4, 
1/5, 2/55, 3/5, 4/5,.3/7;-and.4/7 corresponding’ tom ivy -= 02 1. 25-3) ts 
to be found in a companion report [8]. 

A representative tabulation of the hyperbolic-like LCS functions 


F(x), Hy 0)» and TG) is given in, for example, Tables 8A and 8B of 


the Appendix for a= 3/5. The values of the argument x are the same as those 


used for the tabulation of the hyperbolic functions by Abramowitz and Stegun 


on 


[1]. We observe from Table 8B and (4.13) that the limiting value of 

TGs) as x>+o (here a = 3/5) is quickly reached, with three-decimal- 
place accuracy already attained for x= 4.5. Moreover, the use of these 
tables (specifically, Tables 8A and 8B of the Appendix) for combat analysis 


is illustrated in the next section. 


7. Numerical Examples 


In this section we examine a couple of numerical examples to show 
some of the insights that may be gained into the dynamics of combat between 
two homogeneous forces from our results (see also [6]). These examples 
illustrate the use of the LCS functions Fix), Hy) and T Gx) for 
analyzing "aimed-fire'" combat modelled with the power attrition-rate coefficients 
with "no offset" (5.1). As in [4-7], we consider S. Bonder's model (3.2) 
for the constant-speed attack against a static defensive position. We will 
focus on the use of the LCS functions for predicting force annihilation, 
since the computing of force-level trajectories with Lanchester functions 


is adequately handled elsewhere (see [4-5])., 


oe 


Let us accordingly consider the constant-speed attack of a homogeneous 
Y force against the static defensive position of a homogeneous xX _ force 
(see Section 3 above for further modelling details, especially Figure 1). 
For our numerical computations, we assume that the fire effectiveness of the 


Y weapon system varies linearly with range, i.e. 
a(r) = (7.1) 
0 for :R. <=... 
eo 
and that the fire effectiveness of the X weapon system varies quadratically 


with range, i.e. 


B(r) = CJa2> 


with RY =R i.e. both weapon systems have the same maximum effective range. 
In other words, p = 1 in (3.4) and vy = 2 for &8(r). We consider a battle 
modelled by the input data given in Table II. In terms of time as the 
independent variable, the attrition-rate coefficients (7.1) and (7.2) become 


vie (3.3) 


a(t) = k(t + C) and b(t) = k(t + ey". (7.3) 
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Table II. Input Data for Numerical Examples 


0.06 X casualties/minute/Y firer 


eR 
i 


0.6 Y casualties/minute/X firer 


R = R. = 2000 meters 
a B 


5 miles/hour 


< 
i] 
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where Ry= Rg, 


Re 


eee ae ye , and ek -«(Z) (7.4) 
From the input data given in Table II, we compute the parameter values shown 
in Table III, since the transformed X force-level equation is given by (3.9) 
with q = (v+1)/(V.+v+ 2), p = 1-q, » = 1, and v = 2. Thus, the X force 
level may be computed with FC) and HD) with a =q = 3/5. Force- 
annihilation prediction involves the limiting value of T(t) = H,_ (I/F) 
as t>+”, From Table 8B of the Appendix and Table III, we note the predicted 


agreement between P(1-0)/T(a) and the limiting value of T, (x) a6 3 or 


[recall (4.13)] for oa = q = 3/5. We now consider two cases: (I) R, = 2000 


0 
meters, and (IT) Ry = 1250 meters. 
When Ry = 2000 meters (see Figure 3 of [4]), we have C = 0 and 
cae 0. The maximum time that the battle can last is tax = Ro/v = 14.91 minutes, 


since at this time the attackers reach their final objective, i.e. the 
defender's position (again, see Figure 1). We now consider the qualitative 
behavior of the w=1, v= 2 force-level trajectory shown in Figure 3 of 
[4]. Theorem 2 tells us that the X force can be annihilated in finite time 


if and only if 


Xo r(p) dy q-P 
ee BEES, 
Yg FP Cq) Op = , (7.3) 


where q = 3/5 and p= 1-q. Using the numerical values in Table III, we 
compute from (7.3) that the X force can be annihilated in finite time if 
and only if 
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Table III. Parameter Values for Numerical Examples 


kK = 4.0233 
k, = 2.6979 
p= 2/524 
T'(p)/T(q) = 
A= 0 


x 


1 


1073 X casualties/(minute)"/y firer 


10> Y casualties/ (minute) ’/X firer 


3/5 


48951 
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meee ATEN 4 (7.4) 


When the X _ force can be annihilated, its annihilation time is given by 


(5.8), which we rewrite here as 


h p-q 
) es! I 
E Aete y= —- — re ae , (7.5) 
"AR 


where 
2x 
( I _ (utvt2) /2 (7.6) 


Thus, for the numerical values given in Table III, the time of annihilation 


of the X force is given by 


x 

7 (r(e)) = 3.544 2 , (7.7) 

q a Yo 
with q = 3/5. We will now illustrate further computations for Xo = 10 and 
Yo = 30. From (7.4) we see that the xX force can be annihilated in finite 
time (but we must verify that . Gag) * In this case (7.7) becomes 

Xx 
cee) = 1.18122. (7.8) 


We must now determine t(t5) such that t(t?) = T (1.18122) by using inter- 


polation methods and Tables 8A and 8B. From Table 8A, we find 


ea | 


ae = 1.18172 for = 1.01 
re = 1.17630 for = 1.00 
so that using linear interpolation, we obtain 
xX 
t(t,) = 1.009 , (7.9) 
whence use of (7.6) yields 
X , 
to = 14.24 minutes, C7.10) 


which is less than Bae = 14.91 minutes so that the defending X force is 
indeed annihilated before the attacking Y force reaches its final objective. 
Since r(t) = Ry - vt, we find that force separation at the instant of 


annihilation of the xX force is 


0-6 NaEaeea. . Gan) 


Further results may be computed in a similar fashion and are given in 


Table IV. 
When Ry = 1250 meters (see Figure 3 of [5]), we have C = 5.5923 
minutes, tT, = 0.0975, and t = R./v = 9.32 minutes. In this case 
0 max 0 


Theorem 2 tells us that the xX force can be annihilated in finite time 


if and only if 
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Table IV. Annihilation of the xX Force as a Function 


of the Initial Force Ratio for Ry = 2000 meters 

(x,/y,) — (minutes) (meters) 
OO a a 

0.333 14.24 89.8 
04250 11.61 443.2 
0.200 10.19 633..2 
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2 mpachy 
0 ( A )* P rq) tFyCo) P(e) #9} (7.12) 
utv+t 2 T'(q) = Ale pP) . 

fF (0) P(q) H(t) } 


with q = 3/5 and p =1-q. Using linear interpolation, we obtain from 


Tables 7A and 8A of the Appendix that for the numerical values of Table III 


a = 1.006 , Ege = 0.044 , 
(7213) 
R(t) ="2 3004), Bo (¢.) ="0223 -; 
q 0 p 0 
so that (7.12) says that the X force can be annihilated if and only if 
*0 
So: = Oy 362... (7.14) 
Yo 
When the xX _ force can be annihilated, its annihilation time is given by 
(5.4), which we rewrite here as 
—————— ——————————— + 
7 (- Py 22 rt Ba? 
fae 
T (r(t*)) = - (7.15) 
q a Xo ri P-q 
Et.) 2) o—- H(t.) 
qo Ya utv+t2 q 0 
Yo'AR 
whence for the data of Table III 
x 3.565u, PO2223 
Prt = (7.16) 


0.156u, + 1.004 ° 


where Uy = Xo/Yo- Let us also record here that (3.11) yields 
2/ (utvt2) 
+ vt 2 
t = (espe) aS (7.17) 
I 
We will again illustrate further computations for Xo = 10 and Yo = 30. 


From (7.14) we see that the X force can be annihilated in finite time (but 


j ; : X . 
again we must investigate whether or not e., oe ee In this case (7.16) 


becomes 
Xx 
utr) = 1.533651 ,. (7.18) 


whence Table 8A of the Appendix and linear interpolation yield 


X 
r(t,) = 1.397 , (7.19) 
so-that by (7.17) 
X . 
= 10.63 minutes . (7x20) 
Since t = R,/v = 9.32 minutes < t, we see that the attacking Y force 
max 0 a 


overruns the defender's position before annihilation of the xX force occurs. 


Thus, the battle ends with x, = x(t_) > 0 and Ve > 0 at t.=¢f 


9.32 minutes. Corresponding to te = 9.32 minutes is = 1.1318, and 


then Table 8A of the Appendix yields 
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F(t, = 1.1318) = 1.589 , ave nor oe = 1.973), C7521) 
whence via (2.4), (4.8), (4.9), and (7.13) we obtain 
x. = x(t_) = x(r = 0) = 1.35. (74,22) 


Some further numerical results are given in Table V. Again, these parametric 
results should be contrasted with the single u=1, v= 2 _ force-level 


trajectory shown in Figure 3 of [5]. 


8. Final Remarks 

In the previous section above, we have seen how the LCS functions 
allow one to conveniently obtain much valuable information about the model 
(2.1) with power attrition-rate coefficients (3.8) without having to explicitly 
compute the entire force-level trajectories. Previously we were limited to 
computing only force-level trajectories (see [4-5]). With the availability 
of these tabulations of LCS functions (see the Appendix of this report and 
[8])}, we can now tell who is going to be annihilated and when this event will 
happen without having to compute the trajectories. Not only did we answer 
questions about the qualitative behavior of the model (e.g. force annihilation) 
for specific values of, for example, initial force levels but also for a 
range of values of the initial force ratio (i.e. parametric analysis of 


model behavior). 
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Table V. Annihilation of the xX Force as a Function 


of the Initial Force Ratio for Ry = 1250 meters 
(x./y,) (minutes) ra (meters) 
Xo! VQ a 7 a 
0.333 10.63 t 
0.250 7.56 235.9 
0.200 6.17 422.8 
Tr = 9.32 minutes and «x. = x(r=0) = 1.35. 
max £ 
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The results of this report may be used for other parametric 
analyses, e.g. parametric dependence of battle outcome on attrition-rate 
coefficients. Thus, the contents of this report allow one to develop 
important insights into the dynamics of combat between two homogeneous 
forces with temporal variations in fire effectiveness. With the avail- 
ability of tabulations of the LCS functions, one can now analyze such 
combat modelled by the power attrition-rate coefficients (3.8) with some- 
what the same facility as he can for the constant-coefficient case and thus 
aid in parametric analyses. For further discussions of the significance 
of such results for military operations research, the reader is directed 


to [6] and [7]. 
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APPENDIX: Tabulation of the LCS Functions F (x), H,_\(x), and T(x) for 


l-a 
w= 1/2. 1/3, 2/35, 4/4. Slee 1 by 2/563) oe a) oe 2) 
4/7, S115 479, S195 3/21 yg 31s 6/1 eyo te es 


S/N] 5 12/7. Shel, dnd 16/21, 
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